Off-shell pairing correlations from meson-exchange theory of nuclear forces 
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We develop a model of oflt-mass-shell pairing correlations in nuclear systems, which is based on 
the meson-exchange picture of nuclear interactions. The temporal retardations in the model are 
generated by the Fock exchange diagrams. The kernel of the complex gap equation for baryons is 
related to the in-medium spectral function of mesons, which is evaluated nonperturbatively in the 
random phase approximation. The model is applied to the low-density neutron matter in neutron 
star crusts by separating the interaction into a long-range one-pion-exchange component and a 
short-range component parametrized in terms of Landau Fermi liquid parameters. The resulting 
Eliashberg-type coupled nonlinear integral equations are solved by an iterative procedure. We find 
that the self-energies extend to off-shell energies of the order of several tens of MeV. At low energies 
the damping of the neutron pair correlations due to the coupling to the pionic modes is small, but 
becomes increasingly important as the energy is increased. We discuss an improved quasiclassical 
approximation under which the numerical solutions are obtained. 



I. INTRODUCTION 

The lowest order in the interaction mean-field treat- 
ments of the nuclear pairing in bulk nuclear matter and 
in finite nuclei are conceptually simple and have been rea- 
sonably consistent with the phenomenology of these sys- 
tems. The Hartree-Fock-Bogolyubov (HFB) theory with 
zero- or finite-range (e.g., the Skyrme and the Gogny) 
phenomenological forces has been the standard tool in 
the studies of the nuclear structure to describe the bulk of 
the experimental data [J. At the same time, the calcula- 
tions of the pairing gaps in bulk nuclear matter based on 
the Bardeen-Cooper-SchriefFer (BCS) theory with inter- 
actions modeled by the bare (realistic) nucleon-nucleon 
forces and quasiparticle spectra in the effective mass ap- 
proximation are not inconsistent with the cooling simu- 
lations of neutron stars or the phenomenology of their 
rotation dynamics. Although there is no direct evidence 
for higher order correlations in the phenomenology of nu- 
clear systems, these cannot be discarded on the basis of 
purely phenomenological arguments. The effective forces 
used in the HFB calculations are fitted to the pairing 
properties of (some) nuclei, consequently the contribu- 
tions beyond mean field are built-in automatically (e.g., 
the effective Skyrme or Gogny forces include the contri- 
butions from the ladder diagrams). The values of the 
pairing gaps deduced from the phenomenology of neu- 
tron stars depend sensitively on the physical input, such 
as the equation of state, the matter composition, and 
its dynamical properties. Provided the uncertainties in 
these ingredients the phenomenological values of the gap 
can be off by an order of magnitude. From the theory 
point of view the higher order correlations might not be 
manifest due to large cancellation between the vertex and 
propagator renormalization. 

The on-mass-shell properties of the pairing gaps are 
well understood at the mean-field level. Upon ignoring 
the renormalization effects, the on-mass-shell gaps can 
be directly computed from the nucleon-nucleon scatter- 
ing phase shifts. Since the pairing is a low-temperature 



phenomenon, a well-defined Fermi sea implies that the 
on-shell self-energies of weakly excited quasiparticles can 
be expanded around their Fermi momentum. The magni- 
tude of the gap is then controlled by the renormalization 
of density of states at the Fermi surface. Since the ratio 
of the momentum-dependent effective mass to the bare 
mass is less than unity the reduction of the density of 
states causes a reduction of the pairing gap. 

The early work to construct renormalization schemes 
for strongly interacting fermions introduced wave- 
function renormalization factors in the kernel of the gap 
equation and modeled the pairing interaction in terms 
of particle- hole irreducible vertices These were re- 
lated to the phenomenological Landau Fermi liquid pa- 
rameters which can (in principle) be fixed by compari- 
son with experiments. The Landau parameters can be 
obtained from ab initio microscopic calculations as well, 
which are needed, for example, to describe isospin asym- 
metric nuclear systems that are not directly accessible 
in laboratories ^ IE 1^ - Various aspects of the pair- 
ing in the infinite nuclear matter beyond the mean-field 
(BCS) approximation were addressed during the past two 
decades by employing either more advanced schemes to 
calculated the self-energy and vertex corrections and/or 
modern (phase-shift equivalent) interactions. This work 
includes the microscopic calculations of the effective in- 
teractions and Landau parameters using variants of the 
Babu-Brown 0] or related (e.g., polarization potential) 
schemes of summing the polarization graphs in neutron 
matter 0, 0, ^3 ■ The polarization effects were found to 
suppress the on-shell pairing gaps by large factors. Ear- 
lier, similar reduction was found from the correlated basis 
variational calculations based on the Jastrow ansatz for 
fermionic wave functions 1^ 0, ^3 ■ More recent work 
concentrated on the calculations of the wave-function 
renormalizations and their implementations in gap equa- 
tions employing self-energies derived within Brueckner 
theory at first [3| and second order in the hole-line 
expansion. While again a suppression was found due to 
the wave function renormalizations, supplementing these 
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self-energy corrections with appropriate vertex correc- 
tions largely compensates for the latter effect |0. The 
finite temperature T matrix was utilized to deduce the 
pairing correlations that include full spectral information 
on the system by using the Thouless criterion 17], wave- 
function rcnormalizations in the gap equations, or the T 
matrix in the superfluid phase Iq . The reduction of the 
critical temperature of superfluid phase transition found 
in these approaches is due to dressing of propagators 
in the particle-particle channel, i.e., the analog effect of 
the self-energy renormalizations in the Brueckner scheme 
0, 0. (It is unrelated to the reduction mentioned 
above caused by the polarization of the interaction that 
follows from the resummations in the particle-hole chan- 
nels.) Renormalization group approach has been used to 
resume the particle-hole channels and deduce the pairing 
gap from the matrix elements defined on the Fermi sur- 
face 01; the overall magnitude of the reduction of the 
gap is consistent with the results mentioned above. The 
pairing gaps computed from the renormalization group 
motivated interactions are in good agreement with the 
calculations based on the original bare interactions |23| 
or the effective Gogny interactions 0| below the renor- 
malization scale. 

The off-mass-shell physics of the nuclear pairing is 
much less understood. Contributions to the (complex) 
pairing gap can be generated by time nonlocal inter- 
actions. The simplest example is the fermionic Fock- 
exchange diagram in a boson-exchange interaction model. 
The generic theory of the boson-exchange superconduc- 
tors is known in the form first developed by Eliashberg 
[2_^ to describe the superconductivity in materials where 
the electron-phonon coupling is of the order of unity (see 
also the general references |23,|3l2^)- Because the the- 
ory contains another small parameter - the ratio of the 
phonon Debye energy to the Fermi energy of electrons - 
the vertex corrections (according to the Migdal theorem) 
can be neglected. For the same reason, the quasiclassical 
approximation can be used where the on-shell energy of 
the original Green's functions are integrated out. In our 
context the coupling of the pionic modes to baryons is 
not small and the characteristic energies of the off-shell 
excitations are not small compared to the Fermi energy. 
Below, we shall partially overcome these difficulties, first, 
by including a vertex renormalization by the short-range 
correlations in the pion-nuclcon vertex and, second, by 
applying an improved version of the quasiclassical ap- 
proximation. 

Pairing in nuclear systems, induced by retarded inter- 
actions, was discussed previously in Refs. ^6„ 27, 28] . 
In the context of laboratory nuclei the retarded pairing 
interaction is driven by an exchange of virtual phonons 
- low-lying collective bosonic excitation JSJ ■ In the 
astrophysical context of neutron star crusts the attrac- 
tive pairing interaction is mediated by an exchange of 
real phonons of the nuclear Coulomb lattice of the crusts 
[2^. The approach taken here differs from the models 
above in the nature of the bosonic modes that mediate 



the pairing interactions. While the density fluctuations 
are responsible for the pairing interaction in the models 
above, we concentrate here at the modes that have the 
quantum numbers of pions. 

In what follows, we construct a model of pairing in 
neutron matter that is based on the meson-exchange pic- 
ture of the nuclear interactions. Specific to our model 
is the explicit treatment of the mesonic degrees of free- 
dom including retardation. Unlike the potential models, 
the interaction incorporates medium modifications of the 
dispersion relation of the mesons (pions) and they are 
treated dynamically, i.e., only in the static limit the re- 
tarded interaction reduces to a meson (pion)-exchange 
potential. We note that the meson-exchange picture of 
pairing is known from applications of the Walecka mean- 
field model to the pairing problem 29, .^H, ■ These 
approaches, however, do not include the pion dynam- 
ics or retardation. In computing the renormalization of 
the pionic propagator and the pion-nucleon vertex, we 
shall use a simplified version of the resummation of the 
particle-hole diagrams by modeling the residual interac- 
tion in terms of a contact approximation to the g' Lan- 
dau parameter. This corresponds to the random phase 
approximation (RPA) renormalizations of pion propaga- 
tors in Migdal's theory of finite Fermi systems [33 ■ The 
retardation effects are included by evaluating the Fock 
(exchange) self-energies of neutrons which are coupled 
to dynamical (off-mass-shell) pions. We determine nu- 
merically the resulting off-mass-shell gap function and 
wave-function renormalization sclf-consistcntly by solv- 
ing the coupled integral equations for these two complex 
functions. 

To elucidate our model, wc start with the TTQ-exchange 
interaction among neutrons in a low -density neutron 
matter relevant for the physics of neutron star crusts. 
In the density regime (0.002 — 0.5) poj where po = 0.17 
fm~'^ is the nuclear saturation density, the matter in neu- 
tron stars is an admixture of unbound neutron liquid and 
a Coulomb lattice of nuclei with a c harg e- neutralizing 
background of relativistic electrons [33, llj, 123 . We shall 
assume that the neutrons form a homogenous Fermi liq- 
uid which is characterized by an isotropic Fermi surface. 

The remainder of the paper is organized as follows. 
Section. II is devoted to the formulation of the model in 
the framework of the real time Green's functions formal- 
ism. Section. Ill describes a method of solution of the 
coupled, nonlinear integral Eliashberg equations. The re- 
sults arc discussed for the case under study - the pairing 
in the low-density neutron matter - where the retarded 
pairing interaction is mediated by an exchange of neu- 
tral pions. Section. IV contains our conclusions. In the 
Appendix we discuss the quasiclassical approximation to 
the propagators and its improved version, which is used 
in our numerical work. 
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II. MODEL 

We shall formulate the model in terms of the real-time 
contour-ordered Green's functions (3^ . The imaginary 
time Matsubara formalism is well suited for our purposes, 
since we assume the system to be in equilibrium, however, 
the numerical solution of the equations formulated on 



the real frequency axis are technically simpler compared 
to the solution of the imaginary time equations. The 
latter require a summation of finite number of discrete 
Matsubara frequencies with a certain cutoff frequency 
which is followed by an analytical continuation to the 
real axis; the former are solved directly on the real axis. 



The real-time treatment of a superfluid system requires a 4 x 4 matrix of propagators and self-energies in general 
[3^. A 2 X 2 structure is needed to account for the anomalous correlations due to the pairing. The elements of the 
Nambu-Gor'kov matrix propagator G^p{x, x') are defined as 



-FIp^X^x') Gal3{x,x') 



^i{T^^{x)^l{x')) 
{iUx)i:Ux')) 



{lpa{x)lp0{x')) 



(1) 



where ipaix) are the baryon field operators, x is the space-time coordinate, a and (3 stand for the internal (discrete) 
degrees of freedom, T and T are the symbols for time ordering and inverse time ordering of operators, respectively. 
Each of the propagators above can be viewed as a 2 x 2 matrix in the Keldysh space, if we map the time structure of 
the propagators on a real-time contour. For example, the elements of the normal propagator in the Keldysh space, 
Gap{x, x') (the hat hereafter indicates that the propagator is contour ordered) are defined as 



G^f,{x,x') G^p{x,x') 



-i{M=c)^^lix')) -i{f^c.{x)^l{x')) 



(2) 



The off-diagonal elements of the matrix |(2J) have their time arguments on the opposite branches of the time contour, 
while the upper (lower) diagonal elements have their both time arguments on the positive (negative) branch of the 
contour (hence the time-ordering symbols; see for details Ref. |3y|)- The remainder propagators are defined in a 
similar way. We note that even in the non-equilibrium theory only a limited subset of the elements of the matrices 
above carry new information. The Keldysh structure will be used below as a convenient tool for writing down the 
finite temperature expressions for diagrams. Eventually, wc will find that the final expressions can be written only in 
terms of the normal and anomalous retarded propagators. 

The 4x4 matrix Green's function satisfies the familiar Dyson equation 
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G^Ax^x') = G^Ax^x') + > ^ /d^x"d4x"'G„ ra:,a:"')S .(x"',x")G,,(x",a;') 



(3) 



where the free-propagator matrices are diagonal in both spaces and the matrix structure of the self-energies is identical 
to that of the propagators (see for details Ref. 0|). Fourier transforming Eq. Q with respect to the relative 
coordinate x = x — x' one obtains the Dyson equation in the momentum representation; in this representation the 
components of the Nambu-Gor'kov matrix obey the following coupled Dyson equations 



Gapip) = Goaisip) + Goa-f{p) T.^s{p)Gsi3{p) + ^-is{p)Fhp{p) 

F^M - G^^-,{^p)[K\,{p)Gsp{p) + ^is{-p)Fhp{p)\. 

Faf3{p) = Goayip) ^^s{p)Gsp{~p) + ^-,&{p)Fsfl{p) , 

Gaf3{-p) = Goapi-p) + GQa^{~p) 1:.ff,{~p)G Sfii^p) + ^\s{p)Fsfj{p) 



(4) 
(5) 
(6) 
(7) 



where p is the four momentum, Gajsip) and Goq/3(p) are the full and free normal propagators, F^ ^/^{p) and Fap{p) are 
the anomalous propagators, and S„^(p) and A(j^(p) are the normal and anomalous self-energies; the greek subscripts 
are the spin/isospin indices; summation over repeated indices is understood. We shall assume that the interactions 
conserve spin and isospin, i.e., Gapip) = SapG{p) and Tjap{p) — Sap'S{p) and concentrate below on the pairing in the 
state of total spin S = 0, total isospin 1=1, and orbital angular momentum L = 0. Thus the anomalous propagators 
and self-energies must be antisymmetric with respect to the spin indices and symmetric in the isospin indices 



F'fapip) = gapF'i{p), Fap{p) = gapF{p), 

^ipiP) = gapA^p), Aap{p) = gapA{p), 



(8) 
(9) 
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where gaf3 = icTy is the spin matrix with ay being the second component of the vector Pauh matrices; the unit matrix 
in the isospin space is suppressed. The Dyson equations Q)-l(7J) can be written in terms of auxihary Green's functions, 
which describe the unpaired state of the system 



GU^p) = GoM-p) + G^^{-p)t^s{-p)Ggp{~p). 
Combining Eqs. lEJ-IZJ) and (|1(J|I and Hll|) we find an equivalent form of the Dyson equations 

= G^.^{p)[s^p + A^s{p)FU0{p) 

FKpip) = G^^i-p)Alsip)Gspip), 
hp{p) = G^^{p)k^s{p)Gsp{-p), 
Go,p{-p) = G^.^{~p)\5-,p + K\,{p)Fsp{p) 



(10) 
(11) 

(12) 

(13) 
(14) 

(15) 



The time- reversal symmetry implies that Aa/3(p) = 
[A^^(p)]*, i.e. the pairing correlations in the system 
are fully specified if wc know one of the two anomalous 
Green's functions; in other words, it is sufficient to solve 
Eq. ((T^ with either Eq. ifT^ or Eq. ifTH) . If we are in- 
terested in the equilibrium properties of the system, it is 
convenient to solve Eqs. (|12|) - (|15|l for the retarded prop- 
agators. The form of these equations remains the same as 
we replace the contour-ordered functions by the retarded 
(more precisely, the contour ordered Dyson equation is 
invariant under unitary transformations which bring the 
Keldysh matrix to the so-called 'triangular' form in which 
the diagonal elements are the retarded and advanced 
propagators ; the simplest example of such a transfor- 



mation is given by the unitary matrix [/ = (1 



■ la.. 



,)/^/2 



which acts in the Keldysh space). We first substitute in 
Eqs. (|10() and (|11() for the normal state retarded propaga- 
tors G'^f{±p) = 5c.p[±{uj + ni)-^p-T.^{±p)]-\ where 
is the free quasiparticle spectrum in the normal state. 
After decomposing the retarded self-energy into its even 
and odd in uj components, S^(p) — Yj^{p) -I- S^(p), we 
further define the wave-function renormalization Z{p) — 
1 — uj^^Y.^{p) and the renormalized quasiparticle spec- 
trum ^* — + ^^(p, 1^*). With these definitions, the 
solution of algebraic equations H12|) and (|14|l is 



G«(p) 



F«(p) 



^Z{p)+C, 



{Lj + i7j)^z{p)^-^;^-AR{py 



(cO + i7?)2Z(p)2-^;2„^fl(p)2' 



(16) 



(17) 



where we used the relation Aq,^A^^ — — A^. For the 
purpose of analytical representation of diagrams for the 
self-energies it is convenient to introduce additional real- 
time Green's functions by the relations 



G><ip) = *^G(p)/>'<(t^), 
F>^<{p) = zAF{p)f>^<{u;), 



(18) 
(19) 



where Ag{uj) and Af{uj) are the normal and anomalous 
spectral functions and /^'^(w) are the Wigncr distribu- 
tion functions. In nonequilibrium theory these propa- 
gators are the off-diagonal elements of the Keldysh ma- 
trix |(2Jl (i.e., have their time arguments on the opposite 
branches of the timecontour) and are the solutions of 
transport equations. However, in equilibrium the rela- 
tions (|18|) and l|19|) do not contain additional informa- 
tion, for the spectral functions Ag(uj) and Af(uj) can 
be related to the retarded propagators by the their spec- 
tral representation (here we drop the three-momentum 
dependence of the functions) 



G«(^) 
F«(c.) 



°° de Agje) 
°° de Apie) 



2tt uj — e - 



irj 



(20) 
(21) 



and the Wigner distribution functions reduce to the 
Fermi distribution functions for particles and holes: 
f<{u) = f{Lo) = [expiPu) + l]"i and f>{u;) = f{Lo) - 1, 
where (3 is the inverse temperature. In other words, the 
relations H18|l and H19|l establish a one to one correspon- 
dence between the retarded and G^'^(p), F^''^{p) prop- 
agators. 

To model the strong interaction, we separate its long- 
range and short-range components. The long-range part 
of nuclear interaction is dominated by the pion dynamics 
with the characteristic scale set by the pion Compton 
wavelength = 1.4 fm. The nonrelativistic form of 
the effective pion-nucleon interaction Hamiltonian is 



(22) 



where </> is the pseudoscalar isovector pion field satisfying 
the Klein-Gordon equation, /tt is the pion-nucleon cou- 
pling constant, is the pion mass, cr and r are the 
vectors of Pauli matrices in the spin and isospin spaces. 
For static pions the one-pion-exchange interaction among 
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neutrons in the momentum space has the form 



J TT 



[ti • 0-2 + ^12 (n)]Ti -Xi, (23) 



where q is the momentum transfer, n = q/q, and 512 (n) 
is the tensor operator. It is known to reproduce the low- 
energy phase shift and, to a large extent, the deuteron 
properties fs?]. Below, the pairing correlations are eval- 
uated from the diagrams which contain dynamical pi- 
ous with full account of the frequency dependence of 
the pion propagators; the static results could be recov- 
ered, and a relation to the phase shifts established, only 
in the limit w — *■ in the pion propagators. The in- 
termediate and short-range dynamics is dominated, re- 
spectively, by the correlated two-pion exchange p and 
other heavy meson-exchanges. Short-range correlations 
are crucial for a realistic description of low-energy phe- 
nomena, since the response functions calculated from 
one-pion-exchange alone would lead to an instability of 
nuclear matter at the nuclear saturation density (pion 
condensation) which is not observed. The short-range 
correlations can be modeled by an effective interaction of 
the form IstII 



Korr(w, g) 



g'{uj,q)(Ti ■ (72 



+h'{u;,q)cri ■ ctj Suin) ti ■ tj, (24) 



where g' and h! are the dimensionless Landau param- 
eters. Equation (|24() contains only the spin-isospin part 
of the short-range interaction, which is relevant for the 
RPA renormalization of the vertices and the polarization 
tensor. The parameters g' and h' describe short-range 
correlations on the scale Zcorr which is much smaller than 
the Compton length of the pion, hence, they should be 
smooth functions of the momentum transfer q and the en- 
ergy of excitations m. Frequently, they are approximated 
as g' = const and h' = (?corr^)^x const. Furthermore, of- 
ten only the contribution from the g' parameter is kept, 
since the contribution of the tensor component ~ h' is 
numerically small. The phenomenological features of the 
pion-nucleon systems described above are consistent with 
the microscopic theories 0, • The full meson propaga- 
tor obeys the Dyson equation 



b{p) - boip) + ba{pmp)b{p), 



(25) 



where Do {q) is the free meson propagator and the polar- 
ization tensor is defined as 



n(g) = -Tr 



d'^p 
(2^ 



ito{q)Gip + q)Gi-p)tiq),i26) 



where T^^q) and T{q) are the bare and full pion-nucleon 
vertices; To{q) = —{fTr/m.^) {cr ■ q). The full pion- 
nucleon vertex is defined by the integral equation 

r(9) = f o(g) + Tr / i f i(g)G(p + q)G{-p)t{q), 



(2^) 



(27) 



where Ti{q) is the short-range part of the particle- hole 
interaction, which can be approximated by Eq. 124|) . In 
general the propagators in Eq. H26|l and 127|l are 2x2 
matrices in the Nambu-Gor'kov space. The contributions 
from the anomalous sector, however, involve propagator 
products which contain powers of F ^ each contributing a 
suppression factor 0(A//i), where /i is the Fermi energy. 
Therefore, we shall keep below only the contribution from 
the nonsuperconducting propagators to the polarization 
tensor, i.e., the renormalization of the pion dispersion re- 
lation is carried out in the normal state. As in the baryon 
sector, we introduce the pion Green's functions D^'^{q), 
which are the off-diagonal elements of the nonequilibrium 
matrix propagator (i.e., have their time arguments on the 
opposite branches of the time contour) by the relations 

D<{q) = ^iB{q)g{u;)~iB{^q)[l+g{~LO% (28) 
D>{q) - -iB{q)[l+g{uj)]~iB{^q)g{^uj), (29) 

where g{Ld) — [exp(/3a;) — 1]"^ is the Bose distribution 
function and B{q) is the pion spectral function defined 
over the frequency range w € [0; oo]. The latter is related 
to the retarded component of the polarization tensor, Eq. 
(|26|l . by the relation 



B{q) 



-2lmn"{q) 



Ren«(q)]2 + [lxnm{q)f 



.(30) 



The one-loop-polarization tensor for neutral pions in a 
noninteracting neutron gas at zero temperature is given 
by 



-<pLi(^,q), 



(31) 



where (j)L{^i q) is the Lindhard function. The RPA renor- 
malization of the polarization tensor then leads to 



n«(c.,q) 



2„2 



1 



2f2 
mi 



4'Lii^,q)- 



(32) 



The form of the polarization tensor thus completely de- 
termines the spectral function of pions, which is the cen- 
tral quantity for the calculations of self-energies. The 
time nonlocal parts of the normal and anomalous self- 
energies are due to the Fock contributions defined as 



±{p) = -Tr 
Mp) - -Tr 



(27r)4 
d-^q 
(2^ 



To{q)D{q)G{p-q)r{q), (33) 
^ro{q)b{q)F{p^q)r{-q).{34) 



Note that only the central part of pion-exchange in- 
teraction contributes to the pairing interaction in the 
5 = 0, / = 1 state. To project out the pairing inter- 
action for these quantum numbers we use the identity 

{(Ti ■ q){(T2 ■ q) = i(cri- 0-2)5^ 



1 



[3{ai ■ q){cT2 ■ q) - (cTi ■ a2)q^] , (35) 
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FIG. 1: Top panel: The baryon Fock-exchange self-energies 
for normal (left graph) and anomalous (right graph) sectors. 
The solid lines correspond to the fermions, the wavy lines 
- to the pions; the blobs schematically indicate the RPA 
renormalized vertices, the dots - the bare pion-nucleon ver- 
tices. Bottom panel: The normal and anomalous Hartree di- 
agrams. The dots stand for the contact Landau interactions; 
the straight lines are shown for clarity. 



drop the second (tensor) term and substitute 
(cTi • (T2){ti ■ T2) — —3. To obtain the fuU 
pion-neutron vertex we use the same procedure as for 
the polarization tensor to include the effects of the 
short-range correlations. However, we do not dress 
the pion-nucleon vertex by pion mediated interactions, 
therefore our model corresponds to the lowest order 
contribution in the number of pion-nucleon vertices. The 
top panel in Fig. 1 shows the Fock diagrams correspond- 
ing to Eqs. Ip!^ and the exchange diagrams, which 
are related to the direct ones by crossing symmetry are 
not shown, but are taken into account in the definition of 
the polarization function (|32() . The bottom panel in Fig. 
1 shows the Hartree diagrams which contribute to the 
mean field and the BCS interaction via the short-range 
component of the interaction (the contribution of the 
pions to these diagrams vanishes). 



J 



The analytical expressions for the retarded components of the self-energies can be read off from the diagrams 



E«(^) = 



A«(c.) 



dLd T,>{uj)-T.<{Ld) 
J 27r to — u + it] ' 
du! dcj' D> {uo')G> [Co - u)') - D< {C))G< {Co - lo') 
-00 27r Jo 27r uj - u) + ir^ 

did A>{lo) - A<iuj) 
J 27r Lu — Lu + irj 

duo' D> {lo')F^> {oo-lo')-D< {Co)F^< (w - uj') 

,2^ 



27r 



Lo — LO -\- irj 



(36) 



(37) 



where we used the dispersion relations for the retarded self-energies S]^(w) and A^(w) [in Eqs. ^ and (jSZl the 
momentum variables and the vertices have been suppressed] . Substituting the spectral representations of the Green's 
functions in the above expression we arrive at 



I]«(^,p) 
A\uo,p) 



Tr 
Tr 



d^q 
(2^ 

d^q 
(2^ 



' de 

'de 
,27r 



To{q)AG{e,p- q)C{LU,s,q)Tiq), 



To{q)AF{e,p- q)C{LU,s,q)Tiq), 



where 



Ciu;,e,q)= I ^Biio',q 



f{e)+9{co') , l-f{e) + g{uj') 



UJ' 



UJ - 



IT] 



(38) 
(39) 

(40) 



The expressions (|38|) - (l40ll for the Fock self-energies are general to the extent that the approximations to the spectral 
function of pions and the pion-nucleon vertices have not been specified. In the numerical computations below we 
shall use the RPA renormalized spectral function of pions and the one-pion-exchange vertices which are dressed by 
the short-range correlations; clearly, our choice of the resummation scheme is not unique, but it is known to include 
the dominant set of graphs that renormalize the modes related to the long-range perturbations. 

Next we proceed to integrate out the dependence of the self-energies on the on-shell energies by implementing 
the improved quasiclassical approximation (IQCA), which is described in the Appendix. Since at low temperatures 
the particle momenta lie close to their Fermi surface, the quasiclassical approximation (QCA), in its common form, 
expands the functions with respect to the small parameter Sp/pp, where 6p = |p — pf| is the deviation of the particle 
momenta from their value on the Fermi surface, and keeps the zeroth order in 6p terms. Thus, the on-mass-shell 
physics is constrained to the Fermi surface because of the degeneracy of the system. Independently, the QCA assumes 
that the typical off-shell excitation energies are much smaller than the Fermi energy, so that the momentum integrals 
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over the Green's functions can be evaluated within infinite hmits. The IQCA avoids the latter approximation since, 
as we shall see below, the ofi-shell excitation energies arc of the same order of magnitude as the Fermi energy. 
Upon implementing the IQCA, Eqs. H38|l and (|39|l can be written as 



(41) 



^'^{PF,L^) ^ - J duj'K{Lu'){g{Lu')[G'^{Lu+uj')+GQ{Lu~Lo')] 
+ p ^Im [G«(£)]J(c.,c.',e)}, 

+ r ^Im[f^Q(£)]J(a;,a;',£)|, 

J-OQ J 

where YP{uj,pf) and /^^{cu^pp) are the quasiclassical counterparts of the retarded self-energies, 



(42) 



Jie. = i^— + ' -/(f . , (43) 

e — Lo — Lo' — ti] e — Lo + Lo' — trj 

and we defined a momentum averaged (real) interaction kernel as 

^(^) = .0^3 I"' dqq ( B{q, Lo) Tr {ro(9)r(q)}. (44) 
{2ttYpf Jo Jo 

Formally, Eqs. (|41|l and 1)42(1 are a coupled set of nonlinear integral equations for the complex pairing amplitude and 
the normal self-energy (i.e. are equivalent to four integral equations for four real quantities). This becomes explicit if 
we note that the Green's functions on the right-hand side of Eqs. (|41|l and (|42() can be expressed in terms of certain 
integrals 

G'^iiv) = ujZ'^{u;)h{Lj) + hiij), F^iiv) = -/S9h{uj), (45) 

which are defined by Eqs. and Ht)4|l of the Appendix and are functions of the pairing field /Sf^ {p f , ljj) and the 

wave- function renormalization Z^{pf,oj). 



Let us briefly comment on the physical content of Eqs. 
(|41|l and l|42|l . They describe the pairing and wave- 
function renormalization among degenerate fermions 
which live on their Fermi surface and interact via a force 
that depends on the momentum transfer. The effective 
pairing interaction K(io) is an average over the momen- 
tum transfers that are accessible to particles constrained 
to their Fermi surface {0 < q < 2pf) with a weight given 
by the pion spectral function (i.e., the probability of find- 
ing a pion of momentum q at a fixed energy lo). Although 
the function K{uj) is independent of momentum, it is 
not equivalent to a contact approximation to the (nonre- 
tarded) S-wave interaction in the potential models; the 
latter would correspond, in the present context, to the 
limit of zero momentum transfer. Let us recall again 
that the interaction kernel, via the momentum averaged 
spectral function of pions, incorporates the modification 
of the pionic spectrum by the nuclear medium - a feature 
that is missing in the potential models. 

We turn now to the anomalous Hartree contribution, 
shown by the second diagram in the lower panel of Fig. 
1 (the first diagram renormalizes the on-shell particle 
mass) . Approximating the 5- wave interaction in terms of 
the Landau parametrization, Fo{uJ,q) + Go{ij-',q){cri-(T2), 
one finds for the mean-field BCS contribution to the pair- 



ing in the spin-zero state 

Aic^(a;,p) ^ -tj -0^[Fo{L.',p')-3Goiu;\p')] 

xF''<{uj + uj',p + p'). (46) 

If the Landau parameters are further approximated 
by constants, the quasiclassical approximation to the 
Hartree self-energy becomes 

/OO T 
— /(u;)Im [fQ{cu)],{A7) 
-OO 

where v{ii) — pFiri* /{2tt'^) is the density of states at the 
Fermi surface. 

In the zero-temperature limit, Eqs. (|41|1 and 14211 take 
the simple form 

Z^iuj) = 1-- / — ImrG^(e)l L(uj,e), (48) 

^ J-oo 

A'^iu) = / —lm[F'^{e)]L{Lu,e), (49) 

J — OO ^ 
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where the effective (complex) retarded interaction is defined as 



duj'K{uj') 



u + u' ~ irj 



IT] 



(50) 



where 6{e) is the Heaviside's step function. Note that the thermal occupation probability of the pions vanishes in the 
zero-temperature limit and they contribute only to the effective interaction L{uj,e). Note also that due to the time- 
reversal symmetry of the problem the energy integrations in Eqs. (|48|l and (|49|) can be restricted to the positive energy 
domain upon using the property A^{—e) = A*^(e), the time- reversal properties of the quasiclassical propagators (see 
the Appendix) and the symmetries of the effective interaction L{uj,e). The form of Eq. (|48|l remains the same, but 
the effective interaction is replaced according to L{uj,e) L_{uj,e), while Eq. H49|l becomes 



AQ(^) ^ {lmk«,(e)l L+iu;,s)+ Im [F<!5^„(e)] L_(c.,e)} 



(51) 



where the subscripts on the propagators refer to their odd and even in the energy variable parts and the new kernels 
are defined as 



L±{uj,e) 



dw'Kiuj') 



1 



1 



+ UJ + Uj' + IT] e — UJ + Uj' — 17] 

The zero-temperature expression for the BCS contribution to the pairing gap becomes 



(52) 



Ages = (^0 - 3GoM^^) r - {im [F^^,^{e)] - Im [F„l(e)]} 



(53) 



The net pairing gaps in the quasiparticle spectrum is the sum of the Fock-exchange and Hartree terms A(Li;) 
AQ(w) + A^pg, defined by Eqs. ^ and 



III. THE ALGORITHM AND RESULTS 

Equations H48|) and H49|l form a coupled set of four non- 
linear integral equations for four real functions: these are 
the real and imaginary parts of the gap function A\{lo) 
and A2(a;), and the real and imaginary parts of the wave- 
function renormalization Z\{lS) and Zi{io). Since the in- 
teraction kernel L{u;,e) is constructed from the spectral 
function of pions, renormalized by the nuclear medium 
in the normal state, it serves as an input to Eqs. H48|) 
and H49|) and can be computed prior to their solution. 
Thus, as a first step, the L{u},e) function is generated 
on a two-dimensional energy grid. The starting point is 
the construction of the pion-spectral function from the 
real and imaginary parts of the RPA renormalized po- 
larization function. The RPA polarization function is in 
turn constructed from the free-particle polarization func- 
tion (or the Lindhard function) which has an analytical 
form both at zero and finite temperatures. From the pion 
spectral function and the RPA renormalized vertex T(q) 
we then construct the momentum-averaged effective in- 
teraction K{uj) according to Eq. (|44|l in a typical energy 
range < w < 200 MeV. The first part of the computa- 
tion is accomplished by constructing the real and imagi- 
nary parts of the function L(u;,e) on a two-dimensional 
energy grid. Typical computations were carried out on 
a 100 X 100 point mesh. The range of energy integra- 
tions was extended up to 400 MeV to ensure that the 



integrands in Eqs. H48|l and (|49|l have dropped to values 
below the accuracy of the integration. Note that the eval- 
uation of Eq. H50|l requires principle value integrations 
due to the simple poles in the integrand and, in addition, 
evaluation of integrals which are singular at the bound- 
ary. 

In the second part of the computation the four inte- 
gral equations represented by Eqs. (|48|l and H49|) are 
solved self-consistently using as an input the values of 
the function L{Lj,e) stored on a two-dimensional energy 
grid. An iterative procedure is used. The iteration starts 
by assigning to the functions the values Ai(u;) = const, 
A2(w) — 0, Zi{ijj) = 1 and Z2{uj) = 0, where the con- 
stant is set equal to the typical scale ~ 1 MeV. The 
newly computed functions are reinserted in the kernels 
on the right-hand side of the equations and the itera- 
tions are repeated until convergence is achieved. Note 
that the kernels of Eqs. (|48|l and (|49|) are singular; in 
the first of these equations the singularity is always at 
the lower boundary of the integration region, while in 
the second equation singularities occur both at the lower 
end point and within the integration regions (the lat- 
ter is understood in the Cauchy sense). The positions 
of the singularities, of course, change from one iteration 
to another. An efficient use of the quadratures can be 
achieved by dividing the integration region into subinter- 
vals: a 24 point Clenshaw-Curtis integration was used on 
the intervals which contain singularities; in the remain- 
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FIG. 2: The frequency dependence of the free (light lines) 
and RPA renormalized (heavy lines) pion spectral function at 
several fixed values of the momentum transfer q/pp ~ 0.5, 
1.0, 1.5 and 2.0. The density corresponds to pp ~ 0.4 fm~^ 
in the left panel and to pp — 0.6 fm~^ in the right panel. 



FIG. 3: The frequency dependence of the averaged over 
the momentum transfer (dimensionless) kernel K{lu) derived 
from the RPA renormalized pion spectral function for pf = 
0.4 (solid), 0.5 (long-dashed), 0.55 (dashed) and 0.6 (dashed- 
dotted lines). 



der intervals the standard Gaussian rule was used. A 
few iterations were sufficient to achieve a convergence to 
accuracy 10~^ on all points of the energy mesh. 

Below, we present the results obtained for zero- 
temperature low-density neutron matter relevant to 
the physics of neutron star crusts. The density is 
parametrized in terms of the Fermi momentum p — 
p|n/(37r^) and we consider the range 0.4 < pf ^ 0.6 fm 
~^ which is below the maximum of the on-shell pairing 
gap at about pi? ~ 1 fm~^. The gap decreases at larger 
densities due to the repulsive component of the pairing 
force and we need to include (at least) p-meson-exchange 
to reproduce this feature. We supplement our input pa- 
rameters with a set of dimensionless (i.e., normalized by 
the density of states) Landau parameters. We keep them 
constant throughout the range of Fermi momenta con- 
sidered. The effective Landau mass is set to m* /m — 0.8 
and we assume g' = 0.6, although we discuss its varia- 
tions towards the end of this section. To evaluate the 
BCS contribution (|47|l we choose the parameter values 
/o = —0.3 and go — 0.3 which are taken from Ref. @. 
These values of Landau parameters satisfy the well-know 
stability criteria. 

Figure 2 displays the frequency dependence of the spec- 
tral function of pions at several fixed momentum trans- 
fers for densities corresponding to pp — 0.4 (left panel) 
and 0.6 fm~^ (right panel). It can be seen that the contri- 
bution of processes with momentum transfers q < pf/'2 
to the spectral function is unimportant compared to the 
contribution from large momentum transfers q > pf- At 



fixed q the spectral function has a Lorentzian shape as 
a function of energy transfer. For very large momentum 
transfers q < 2pF, as the particle- hole excitations with 
higher energies become accessible, the spectral function 
is broadened and its maximum is shifted towards larger 
energies. Compared to the free case, the maximum of the 
spectral functions computed within the RPA is shifted to 
lower energies while the width remains nearly constant. 

Figure 3 displays the frequency dependence of the in- 
teraction kernel K{uj) for several Fermi momenta. The 
shape of the kernel K{uj) is Lorentzian at low and inter- 
mediate energies [characteristic also to spectral function 
B{q,uj) for fixed values of \q\] with an excess contribu- 
tion in the high-energy tail. The contribution from small 
momentum transfers q <C pp to the kernel K{ijj) is neg- 
ligible, since the spectral function for these processes is 
generally small (see Fig. 2) and, in addition, the inte- 
grand in Eq. (|44|l is weighted by the factor q'^, with 
each pion-nucleon vertex contributing a power of q. The 
latter factor cuts off the contribution from small momen- 
tum transfers. The main contribution to the peak of 
K(uj) function originates from the momentum transfers 
\q\ ^ Pf] in this regime the maximum of the spectral 
function coincides with that of the K{ljj) function. The 
excess contribution to the high-energy tail of the K(u}) 
function is due to the excitations with large momentum 
transfers and the magnifying effect of the pion-nucleon 
vertices (a factor (x q^). As the density is increased, 
the K{u}) function extends to higher energies, its maxi- 
mum increase by roughly an order of magnitude over the 



10 



p = 0.4 fm" 




FIG. 4: The e-dependence of the efTective interaction kernels 
L+(u,e) (heavy hnes) and L-{ijj,e) (hght hnes) for fixed val- 
ues of a; = 0, 0.5, and 1 fm"^ for pp = 0.4 fm"^. The solid 
and dashed lines refer to the real and imaginary parts. 



real and imaginary parts of the gap function develop a 
complicated structure: in the low-energy domain, below 
the maximum of the function K{uj) at about 10 MeV, 
the real part of the gap function dominates the imag- 
inary part: |Ai(a;)| > |A2(ti;)|. Beyond the maximum 
the opposite is true up to the energies lo ~ 100 MeV. A 
nonvanishing imaginary part A2(w) implies a finite decay 
time for the pairing correlations due to the coupling of the 
neutron quasiparticles to the pionic modes. The imag- 
inary part of the wave- function renormalization Z2{lo) 
implies, likewise, finite lifetime for the normal quasipar- 
ticle excitation. It is worthwhile to note that the latter 
effect is not independent of the pairing correlations, since 
it emerges from a self-consistent solution of the coupled 
dynamics of the normal and anomalous sectors. The en- 
hancement of the gap and the wave-function renormal- 
ization with increasing density can be traced back to the 
enhancement of the attractive interaction kernel K{lj). 
Including the repulsive components of the interaction by 
treating heavier (notably p meson) exchanges would lead 
to a saturation of the attractive interactions at about 
Pp ~ 1 fm~^ and a depression of the gap at larger den- 
sities. 



Fermi momentum range 0.4 < pp < 0.6 fm^^ with a mi- 
nor shift of its position to the higher energies. Note that 
the observed increase in the interaction strength will sat- 
urate at densities corresponding to pp < 1 fm^^ once the 
repulsive interaction via p meson-exchange will become 
important. 

The panels in Fig. 4 show the energy dependence of 
the real and imaginary parts of the interaction kernels 
L±{uj,e) [Eq. (|52l) ] at several fixed values of the integra- 
tion variable e. Separating the real and imaginary parts 
in Eq. H52|l it is easy to see that (i) the real parts ac- 
quire their maximum when e = uj and (ii) the imaginary 
parts vanish for e > w features that are evident in Fig. 
4. Thus, the main contribution to the kernel of the gap 
equation comes from the 'diagonal' part of the interaction 
kernel L±{e, e). 

The frequency dependence of the real and imaginary 
parts of the gap function and the wave-function renor- 
malization is shown in Fig. 5. The gap function has 
been replaced by a new function A'(a;) = Z~^{uj)A{uj) 
(the prime is dropped below), which is equivalent to a 
redefinition of the energy scale us' = ujZ{uj). The a; — > 
limit of these functions is consistent with what is known 
from the on-shell physics. The real part of the gap func- 
tion Ai {oj — 0) is of the order of MeV and increases with 
density, as is the case for the pairing gaps deduced us- 
ing potential models. The real part of the wave-function 
renormalization Zi{u! = 0) is larger than unity, i.e., there 
is an enhancement of the density of states at the Fermi 
surface. The imaginary parts of these functions vanish in 
the limit w ^ (as they should). Off the mass-shell the 
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FIG. 5: Top panel: The frequency dependence of the real 
(heavy lines) and imaginary (light lines) parts of the gap 
function Ai{lu) and A2(cj) for pp = 0.4 (solid) 0.5 (long- 
dashed) 0.55 (dashed) and 0.6 (dashed-dotted lines) fm~^. 
The on-shell values of the pairing gap are 0.1, 0.7, 1.4, and 
3.7 MeV for pp = 0.4, 0.5, 0.55 and 0.6 fm"\ Bottom panel: 
The frequency dependence of the real (left panel) and imagi- 
nary (right panel) parts of the wave-function renormalization 
Zi{u:>) and Z2{uj). The labeling is the same as in the top panel. 
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FIG. 6: The frequency dependence of the effective interaction 
kernel K{u!) (left graph) and the real Ai(a;) and imaginary 
A2(a') parts of the pairing gap (right graph) for three values 
of the Landau parameter g' = 0.4, 0.6 and 0.8 and pp ~ 0.4 
fm~^. The light lines correspond to the QCA approximation 
with g' = 0.6 for pp = 0.4 fm~-^. 



FIG. 7: The frequency dependence of the real Ai(cj) and 
imaginary A2(w) parts of the pairing gap (top panel) and 
the real Z\{lS) and imaginary ^2(0;) parts of wave- function 
renormalization (bottom panels) for pp = 0.4 fm~^. The 
solid and dashed curves correspond to the results obtained 
with the RPA renormalized and free pion spectral function. 



to values g' ~ 0.6 — 0.8 The experimental evi- 

dence from the studies of spin-isospin modes in finite nu- 
clei suggests that g' ~ 0.6. Figure 6 displays the changes 
in the interaction kernel and the pairing gap func- 

tion A([j) as the Landau parameter g' is varied in the 
range 0.4-0.8. The larger g' the stronger is the attractive 
interaction kernel and, hence, the larger is the pairing 
gap. The overall range of variations of these quantities 
is, however, small and the observables are insensitive to 
the variations of g' . As discussed in the Appendix, the 
ordinary QCA requires that the characteristic off-shell 
energies of the excitation be smaller than the Fermi en- 
ergy a;//i ^ 1 - a condition which is not well satisfied 
for excitation energies in the present model. One finds 
typically ^ 1. The IQCA relaxes the <C 1 as- 
sumption of the QCA, but treats the pairing gaps and 
the wave-function renormalizations independent of the 
momentum variable. Figure 6 (right panel) compares 
the results obtained within the QCA and the IQCA. Al- 
though the qualitative picture is the same in both approx- 
imations, the IQCA leads to smaller gap functions (and 
wave-function renormalizations) because the integrations 
over the on-mass-shell energy are limited to a finite shell, 
while the QCA assumes infinite integration limits. Calcu- 
lations which include full momentum dependence of the 



gap function and wave-function renormalization will be 
needed to gauge the remaining error due to the IQCA. 
Finally, Fig. 7 compares the results obtained with the 
free and the RPA renormalized pion spectral function. 
The shapes of the gap and renormalizations functions 
are qualitatively similar in both cases; however, the RPA 
renormalization enhances the values of the gap at small 
frequencies and, therefore, their on-shell values. More 
generally, according to Fig. 7, the RPA renormalization 
enhances (suppresses) the gap and renormalization func- 
tions in the low- frequency (high-frequency) limit. These 
features are consistent with the frequency dependence of 
the free and renormalized spectral functions displayed in 
Fig. 2. 

IV. CONCLUDING REMARKS 

In this work we proposed a model of pairing in nuclear 
systems based on the meson-exchange picture of nuclear 
interactions that includes retardation effects. The model 
is illustrated on a simplest possible example - neutrons 
interacting via exchange of neutral pions; the short-range 
correlations are included in the spirit of Landau-Migdal 
theory of Fermi liquids. Unlike the treatments based 
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on the potential models, the modifications of the meson 
properties in the medium are taken into account by relat- 
ing the pairing gap function and the wave- function renor- 
malization to the spectral functions of mesons. Wc com- 
puted these medium modifications for the pionic modes 
within the non-perturbative particle-hole (RPA) resum- 
mation scheme, which softens their spectrum, but has no 
qualitative effect on the retardation. The model builds in 
the retardation effect via the coupled Fock-exchange self- 
energies for normal and anomalous sectors. The Hartree 
self-energies receive contributions only from the short- 
range part of the interaction, as pions do not contribute. 
Although we obtain the exact self-energies for our sys- 
tem, these are solved numerically with the improved qua- 
siclassical approximation (IQCA). Since it is difficult to 
estimate a priori the error introduced through the IQCA 
(which fixes the momentum transfer at the Fermi momen- 
tum) our numerical results should be taken with caution. 

The numerical calculations above show that the real 
and imaginary parts of the gap function are of the same 
order of magnitude over a wide energy range. At low 
energies, the real part of the gap function dominates the 
imaginary part, i.e., the Fourier transform of the anoma- 
lous pair-correlation function in the time-domain shows 
a damped oscillatory behavior; the picture is opposite 
in the high-energy regime, typically above tens of MeV, 
where the oscillations are damped out within a 'cycle'. 
Physically, the complex pairing gap implies a finite life- 
time of the Cooper pairs due to the emission/ absorption 
of pions. The real part of the wave-function renormal- 
ization (derived self-consistently with the gap equation) 
is larger than unity, i.e., implies an enhancement of the 
density of states; the finite imaginary part describes the 
damping of the normal quasiparticle excitations due to 
the coupling to the pion modes. The frequency depen- 
dence of the pairing gap and the wave-function renormal- 
ization will have implications on the responses of nuclear 
systems to those electromagnetic and weak probes which 
are sensitive to the energy range where these functions 



vary substantially. 

Extrapolating from the analysis of this paper, it is clear 
that the meson-exchange picture of pairing with retarda- 
tion can be extended to treat problems which have been 
difficult to handle within the potential models. Such an 
extension could incorporate the dynamics of nucleon res- 
onances, for example, to clarify the role of delta's in the 
effective retarded interaction or to deduce their pairing 
properties. Another interesting direction is the appli- 
cation of the model to the strange sector to study the 
pairing properties of the strange baryons that are sta- 
ble in the interiors of compact stars (typically S and A 
hyperons, see Ref. [sM Hal); the role of the kaon dy- 
namics on the pairing properties of the strange baryons, 
etc. The model appears also suitable for treating the 
fermionic pair condensation of the baryons on the same 
footing with the Bose condensation of pions or kaons. 
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Appendix 

The purpose of this appendix is to discuss in detail the 
approximations to the full normal and anomalous Green's 
functions that allow us to remove the explicit momentum 
dependence of the wave- function renormalization and the 
gap function. We start by examining the momentum de- 
pendence of Eqs. (|38|I - I|4U|I in the low-temperature limit. 
Since the formal structure of the normal and anomalous 
self-energies is the same, we shall apply our arguments 
only to the normal self-energy, Eq. (|38|l . Below, we as- 
sume isotropic medium with a spherical Fermi surface. 



J 



Using the spectral representations of the Green's functions (|^ and (PT|l . Eq. can be transformed to the 

following form 



E«(c.,p)= I ^Tr{ro(g)r(g)} j^'^B{u',q) 
-g{uj')G"ip-q,L,-Lo')+ / -J(e,^,c.')ImG^(p-q,e) 



(54) 



where the function J(e, w, uj') is defined by Eq. Directing the axis of a spherical system of coordinates along the 

vector p, the integration measure becomes dPq — q^dq dx dcj), where x is the cosine of the angle between the vectors 
p and q. Upon defining k = p — q and replacing the integration over x by an integration over ^j,, Eq. (|54|l becomes 

S^(a;,Pf ) ^ ""^ / qdqTv{To{q)T{q)} -^B{oj' ,q)H{oj,w'), (55) 
{2ttYpf Jo Jo 27r 
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with 

H{uj,uj') = g{uj') / d^k [G"ik,uj + uj') + G'^{k,uj-uj')] 
dp ri+il) 

— J{e,u,u') d£_klmG"{k,e), (56) 

where = /2m* — /i, the integration Umits arc ^± = {pp ± q)^ /2m* — /i, m* is the effective mass, /i is the chemical 
potential. Without loss of generality, the momentum argument of the self-energy 'E,^{oj,p) is taken as the Fermi 
momentum, i.e., the expressions (|55|l and (|56l) are still exact. 

Quasiclassical Green's functions. Now we turn to the -integrations in Eq. (|56|l . Formally, the quasiclassical 
approximation (QCA) amounts to (i) taking infinite limits ^± — s- ±oo and (ii) approximating the function A(^fe,w) 
and Z{S^k, ^) in the propagators by their values on the Fermi surface A(^(pir), uj) and lo). The integrals over 

the normal and anomalous Green's functions are then computed by a contour integration 

ujZ{uj) + 



da 



(a;)2 ~ £,1- A-R{w)2 -|- i sgn(w)?7 



ITT 




-A«(w) 






+ i sgn 


A(cj)sgn(cj 
r 
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The second approximation [^(fc) ■^(fcF)], assumes that the gap function and the wave- function renormalization are 
smooth functions of momentum, therefore, at low temperatures, they can be approximated by their values at the 
Fermi surface. For the range of possible momentum transfers < g < 2pi? for particles on their Fermi surface, the 
functions need to be constant in the interval k € [0 ; Spp]- Such an assumption is not inconsistent with the results 
of the potential models, but an error estimate will depend on the shape of the chosen interaction in the momentum 
space. Note that in our case, the momentum dependence of the functions is controlled by the shape of the spectral 
function which contributes mainly for q > pf- 

The first approximation (infinite integration limits) uses the fact that the large momentum transfers are important 
in the integration in Eq. H55|l . so that the limits |i^±| ^ /i are large compared to values of the poles of the Green's 
functions in Eqs. H57|) and (|58|l . i.e., the condition ^poic ^ |C±I is satisfied. The main contribution to the q integration 
comes from large momentum transfers pf ^ q < ^pF (see Sec. Ill), so that £,± = {pf ± qY /2m* — /i ~ ±/i [note that 
max (^+) = 8^ and min (f_) = — /i]. In addition, the pion-nucleon vertices contribute a factor of q^ to the integrand 
in Eq. H55|) . which cuts off the contributions from small momentum transfers. The poles of the Green's functions 
in Eqs. ifsTjl and (|55|l are defined by the condition |^poie| = \/'^'^Z{ujY — /S.{ujY < u. The QCA is then applicable 
whenever the condition ^poic/M ^ 1 is satisfied. In our problem, the typical energy scale for the off-shell excitations 
is of the order of 20 MeV which is of the same order of magnitude as the Fermi energies; the parameter ^poie/M could 
be of the order of unity, in which case the conventional QCA will break down. 

Improved quasiclassical Green's functions. Here we shall improve on the second approximation of the QCA by 
keeping the integration limits finite. To decouple the on-shell energy integrations over the Green's functions from 
the remainder of the kernel, we approximate the q-dependent integration limits £±{q) by constants £±{pf)- Setting 
the momentum transfer to q = pf is motivated by the shape of the pion-spectral function which, as noted above, 
contributes essentially in the range pf < q <'^Pf- 

The integrations in Eq. (|55|l and the analog equation for the anomalous self-energy require then two types of 
integrals (^ = ^ hereafter) 



d£ 



UJ 



^Z{u;y -e- A«(cj)2 + i sgn(a;)77' 



(59) 



I u}^Z{uj)^~£^-AR{Ljy + isgn{uj)7j' ^^^^ 



Upon defining a shorthand D — uj^Ziuj)^ — A^((jj)2, the first integral can be rewritten as 



D — S/^ +i sgn{uj)ri D — S^'^ + i sgn{uj)r] 



(61) 
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The first term contains a pole and we apply the Dirac identity 1/a; = P/x — inS^x) to separate the real and imaginary 
parts {P stands for the principal value). The second term is regular and we take the limit 77 — > to obtain 



9{D) 



P 



-D+e 



i 7rsgn(ci;) 



m 



These integrals are standard (e.g. 3.3.21 and 3.3.23 of Ref. '4d\) and the final result is 



(62) 



0{D) 



In 



AD) 



-D) 



-arctanh 



(63) 



The step functions in the last term insure that the poles 
of the propagators are counted only if they lie within the 
integration limits; to obtain the final expression we have 
used the relation arctan(ix) — i arctanh(a:). In the limit 
^± ±00 we recover the right-hand sides of Eqs. H57|) 
and (|58|l . The integral I2 is computed as above and, upon 
using 3.3.22 of Ref. we obtain 



D-e 



[9iD) 



-D)]. 



(64) 



Note that this integral vanishes when the limits are infi- 
nite (^-t ±00) or, more generally, when the integration 



is within finite limits which are symmetrical ( | ^_ | = | C+ 1 ) • 
Under the time reversal the integrals Ii and I2 transform 
according to 



Re /i,2(— t^) = Re /i^2(w), Im Ii{—uj) = — Im 



These transformation properties were used to derive Eqs. 
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